Orbital magnetization and its effects in spin-chiral ferromagnetic kagome lattice 
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Recently, Berry phase in the semiclassical dynamical of Bloch electrons has been found to make 
a correction to the phase-space density of states and a general multi-band formula for finite- 
temperature orbital magnetization has been given [Phys. Rev. Lett. 97, 026603 (2006)], where the 
orbital magnetization A4 consists of two parts, i.e., the conventional part M c and the Berry-phase 
correction part Mn. Using this general formula, we theoretically investigate the orbital magne- 
tization and its effects on thermoelectric transport and magnetic susceptibility properties of the 
two-dimensional kagome lattice with spin anisotropics included. The study in this paper is highly 
interesting by the occurrence of nonzero Chern number in the lattice. The spin chirality parameter 
4> (see text) results in profound effects on the orbital magnetization properties. It is found that the 
two parts in orbital magnetization opposite each other. In particular, we show that M c and Mn 
yield the paramagnetic and diamagnetic responses, respectively. It is further shown that the orbital 
magnetization displays fully different behavior in the metallic and insulating regions, which is due 
to the different roles M c and Mn play in these two regions. The anomalous Nernst conductivity is 
also calculated, which displays a peak-valley structure as a function of the electron Fermi energy. 

PACS numbers: 75.30.-m, 73.43.-f, 72.15.Jf 



I. INTRODUCTION 
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In the semiclassical dynamics of Bloch electrons, the 
Liouville's theorem on the conservation of phase-space 
volume is violated by the occurrence of Berry phase, 
which leads to a modification of the phase-space den- 
sity of states^. This modified phase-space density of 
states, D n (r, k) = l+eB-r2„/ft, enters naturally in the 
semiclassical expression for the expectation value of phys- 
ical quantities, and has profound effects on equilibrium 
as well as transport properties. Here B is the ex- 
ternal magnetic field, and fi„(k) is the Berry curva- 
ture (gauge flux) of electronic Bloch states defined by 
fi n (k) = i(V k u n (k)| x |V k u„(k)) with |u„(k)) being the 
periodic part of Bloch wave for the ro-th band. Based 
on this correction in the semiclassical phase-space den- 
sity, Xiao et ali^ have further derived a general expression 
for the single-particle free energy in the presence of the 
external magnetic field as follows (for a general three di- 
mensional system) 



(i) 

Here /j, is the electron chemical potential, (3 = X/UbT, 
and e* f = £ nk — m n (k) • B is the electron band energy 
in the presence of the external magnetic field. m„(k) 
is the crystal orbital magnetic moment and convention- 
ally defined bj*M m„(k) = -£(e/2fi)(Vk«n(k)| x [H(k)~ 
£„(k)]|Vkii„(k)), where H(k) is the crystal Hamiltonian 
acting on |u„(k)). The orbital magnetization (OM) is 
then given by the field derivative at fixed temperature 
and chemical potential, Ai(r) = — (dF/dH)^ T , with the 



M(r) =J2 J d 3 km„(k)/„(r,k) (2) 
n I 

M c + M , 



where / n (r, k) is the local equilibrium Fermi-Dirac distri- 
bution function for n-th band. The first term in Eq. (2) 
is just a statistical sum of the orbital magnetic moments 
of the carriers originating from the self-rotation of the 
carrier wave packet^, thus we call this term the conven- 
tional part M c of the OM. Whereas the second term Mq 
is the Berry phase correction to the OM. This term is 
of topological nature, arising from a bulk consideration. 
Interestingly, it is this Berry-phase term that eventually 
enters the transport current^. At zero temperature and 
magnetic field the general expression (2) reduces trji£ 



d*k m„(k) 



-J7„(k) [no - e nk ] ) , (3) 



where the upper limit means that the integral is over 
states with energies below the Fermi energy fiQ. As men- 
tioned above, the first term in Eq. (3) is the contribution 
from the intrinsic orbital moment, and the second term 
comes from the Berry-phase correction of the density of 
states. 

Among very few studies on Berry phase effect of the 
OM, Xiao et al^ have investigated the anomalous ther- 
moelectric transport in CuCi^Sci-^Brj; system^ by use 
of their derived relation between anomalous Nernst ef- 
fect (ANE) and anomalous Hall effect (AHE). Ceresoli 
et alj£ have given an extensive tight-binding calculation 
of OM for finite and periodic two-dimensional (2D) Hal- 
dane model^ in a wide range of parameters covering cases 
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of Chern number (a dimensionless integer—) C = and 
C =/= 0. Due to its basic importance in understanding the 
magnetism and transport features of the materials, obvi- 
ously, more work are needed in exploiting the properties 
of the OM in various kinds of realistic physical systems. 
Also the theory itself, such as a full quantum mechanical 
derivation of the OM, remains to be further developed. 

In this paper we extend the study of the OM to the 
strongly correlated electronic systems. More specificially, 
we focus our attention to a typical flat-band ferromagnet 
with spin anisotropies on the 2D kagome lattice. This 
attention is partially motivated by the recently estab- 
lished common point that the spin Berry phase plays an 
important role in the quantum transport in spin-orbit 
couple d 10 ! 11 ' 12 or spin-chira l 13 ' 14 ' 15 ' 16 ferromagnetic sys- 
tems; the latter is exampled by pyrochlore compounds 
R2M02O7 (R=Nd, Sm, Gd), in which the spin config- 
uration is noncoplanar and the spin chirality appears. 
As a consequence, the quantum transport of electrons, 
especially the transverse conductivity a xy is expected 
to be affected by the presence of spin chirality. In 
fact, recent transport experiments on ferromagnetic py- 
rochlores have revealed that the AHE increases as the 
temperature T is lowered and approaches to the satu- 
rated valu o 17 ' 18 . This behavior is intrinsically different 
from the conventional theories 1 ^. One explanation to 
this anomalous feature is that the pyrochlore structure 
has geometrical frustration^ - which consists of corner- 
sharing tetrahedrons. Thus the antiferromagnetic and 
even the ferromagnetic^ interaction between nearest- 
neighbor spins are frustrated. Recently, Ohgushi et al^ 2 - 
have first pointed out that the chiral spin state can be 
realized by the introduction of spin anisotropy in an or- 
dered spin system on the 2D kagome lattice, which is the 
cross section of the pyrochlore lattice perpendicular to 
the (1, 1, 1) direction 2 ^,. In this case, it has been shown- 2 * 
that the presence of chiral spin state may induce nonzero 
Chern number, thus resulting in a quantized Hall effect 
in insulating state, which is expected to have important 
implications to AHE experiments in ferromagnetic py- 
rochlores. 

Motivated by the above work22 on quantized Hall ef- 
fect in the 2D kagome lattice, in this paper we turn to 
study the OM properties and effects in this 2D lattice 
system with nonzero Chern number. We show that the 
OM displays different behavior in the metallic and insu- 



lating regions, which is due to the different roles M c and 
Mq play in these two regions. Also, the ANE and the or- 
bital magnetic susceptibility are investigated, which fur- 
ther illustrate the fundamental role brought forth by the 
Berry-phase contribution and nonzero Chern number. 

II. THEORETICAL MODEL AND CHERN 
NUMBER ANALYSIS 

Now we consider the double-exchange ferromagnet on 
the kagome lattice 2 ^ schematically shown in Fig. 1. 
Here the triangle is the one face of the tetrahedron, and 
the easy axis of the spin anisotropy points to the cen- 
ter of each tetrahedron and has an out-of-plane compo- 
nent. In this situation the three local spins on sites A, 
B, and C in Fig. 1 have different directions and the 
spin chirality emerges. The effective Hamiltonian for the 
hopping electrons strongly Hund-coupled to these local- 
ized spins is given by H = Y^NN^ij rtfy witn *i/ = 
t(Xi\Xj) = te laij cos Here the spin wave function \x%) 

is explicitly given by \xi) — [cos^, e % '^ i sin^j-] T , where 
the polar coordinates are pinned by the local spins, i.e., 
(Xi|Sj|xi) = \ (smtii cos fa, smtfjSin^i, costfj). % is 
the angle between the two spins Sj and S^. The phase 
factor aij can be regarded as the gauge vector poten- 
tial a M (r), and the corresponding gauge flux is related 
to scalar spin chirality Xijk=Si-(S j xSk)— . In periodic 
crystal lattices, the nonvanishing of the gauge flux re- 
lies on the multiband structure with each band being 
characterized by a Chern number 9,24 '. Here the Chern 
number appears as a result of the spin chirality in ferro- 
magnets. Following Reb^ 2 - we set the flux originated from 
the spin chirality per triangle (see Fig. 1) as <f>, which sat- 
isfies = e i ( QAB+aBC+acA ). The flux penetrating one 
hexagon in Fig. 1 is determined as —2(f). We take the 
gauge, in which the phase of t^y is the same for all the 
nearest-neighbor pairs with the direction shown by the 
arrows in Fig. 1. It should be pointed out that the net 
flux through a unit cell vanishes due to the cancelation 
of the contribution of the two triangles and a hexagon. 
Also noted is that the time-reversal symmetry is broken 
except for cases of 0=O,7r. 

The momentum-transformation of the above tight- 
binding Hamiltonian is given by 22 





2cos(k • a^e 1 ^ 3 
2cos(k • a 3 )e- ,; ^ 3 



2cos(k ■ ai)e- ,;< ^ 3 


2cos(k-a 2 )e^ 3 



2cos(k • a 3 )e^/ 3 
2cos(k-a 2 )e- ,; ^ 3 




(4) 



r 



where ai = (— 1/2, — \/3/2), a2=(l, 0), and a 3 = from A to B site, from B to C site, and from C to A site 
(-1/2, V3/2) represent the displacements in a unit cell respectively. In this notation, the Brillouin zone (BZ) is a 
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FIG. 1: (Color online). Two dimensional spin-chiral ferro- 
magnetic kagome lattice. The dashed line represents the 
Wigner-Seitz unit cell, which contains three independent sites 
(A, B, C). It is assumed that each site has a different spin 
anisotropy axis. The arrows means the sign of the phase of 
the transfer integral tij. 



J 



hexagon with the corners of k = ±(27r/3)ai, ±(27r/3)a2, 
±(27r/3)a3, two of which are independent. Note that in 
writing the Hamiltonian (4) we have chosen the unit of 
t cos(i?y-/2) = 1, and set the length of each bond as unity. 
The eigenenergies of the Hamiltonian (4) are given by 



, b(k) 0(k)+27T 
Elk = 4W4r Z COS- 



. b(k) 0(k)-27T 



(5) 



. b(k) 0(k) 



where 6(k)=cos 2 (k 
6>(k)=arg 



ai)+cos 2 (k 
^c(k)+i 



1 



a 2 )+cos 2 (k • a 3 ), 
, and c(k)=cos(k- 



27c 2 (k) 
6 a (k) 

i. Because the relation 



ai)cos(k ■ a 2 )cos(k • a 3 )cos<? 
£ik < ^2k < £3k is always satisfied, henceforth we will call 
these energy bands 'lower', 'middle', and 'upper' bands 
respectively. The corresponding eigenvectors are 



i[4k-W(k.a 2 )] 

\Unk) = G n (k) e l 3 [e nk cos(k • ai) + 2 cos(k ■ a 2 ) cos(k • a 3 )e~^] 
e~ l a [e n k cos(k • a 3 ) + 2 cos(k • ai) cos(k • a 2 )e~^] 

I 



(6) 



where the normalized factors are given by 



G„(k) = 



1 



>/26(k)4 k + (46(k) - 3e 2 k ) cos 2 (k • a 2 ) + 12c(k)e„ k 
I 



(7) 



First let us see the Chern number and Hall conduc- 
tivity of this system, which has been briefly reported in 
Refi^ 2 -. It is clear that the Hall conductivity a xy is equal 
to zero in the time-reversal symmetric cases cf> — 0, n. 
Therefore, we focus on the case of <f> ^ 0, it. In this sit- 
uation there is an energy gap between the two nearest- 
neighbor bands. We assume that the Fermi energy is 
lying in the gap. Then the Hall conductivity is a sum 
over occupied Bloch bands, a xy =(e 2 /h) Y^ CU Cn, where 



the n-th band Chern number is denned by 



c n = ~ [ d 2 kn n {k) 

27r Jbz 



(8) 



= ~ f d 2 kz-(V k x A„(k)) : 
27T Jbz 



where A„(k)=z(it n (k)|VkW n (k)) is the Berry phase con- 
nection (vector potential) for the n-th band. According 
to the expressions for |it„k), one obtains the expression 
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for A„(k) as follows 



A„(k) = 2e„ k G 2 l (k) sin cos fc x 



sm „x 



%/3sin k x y 

To proceed with Eqs. (8)-(9), one may first transform 
the integral of the curl of the vector potential A„ over 
BZ to the line intergal of A„ along the BZ boundary by 
use of Stokes' theorem, and then apply the complex con- 
tour integration technique and residue theorem to sinu- 
soidal functions. After a straightforward derivation, one 
obtains C\ = — sgn(sin0), C 2 = 0, and C 3 =sgn(sin0), 
respectively, which means that the quantum Hall effect 
is realized. However, this purely mathematical calcula- 
tion of Chern number is not favored by theoretical physi- 
cists, who instead resort to the physical connotation that 
the vector potential A n and gauge flux tt n are endowed 



with. Correspondingly, we start this gauge-field analysis 
with the notation that the value of Chern number is in- 
variant under gauge transformation \u' nis )=e l9 ' n ^\u n \ s ) , 
A^(k) = A„(k) — V k g„(k), where g n (k) is an arbitrary 
smooth function of k. The gauge of the wave functions 
Eq. (6) has been chosen in such way that the first com- 
ponent is real. If this gauge choice is applicable in the 
whole region of the BZ, the Chern number will obviously 
be zero. However, at two equivalent BZ edge points 
k = (0, ±7r/-\/3) on BZ boundary, one can find that 
the wave functions |u„ k ) m Eq. (6) is ill-defined since 
both its denominator and numerator are zero at these 
two points. This means that the used gauge cannot ap- 
ply to the whole BZ and one needs to render a gauge 
transformation to avoid the singularity at ko. For this 
one writes down the other set of eigenvectors as follows 



J 



e [e„k cos(k • ai) + 2 cos(k • a 2 ) cos(k • a 3 )e 1, 
(V;,!k) j i[4 k _4cos 2 (k.a 3 )] 

[e„ k cos(k • a 2 ) + 2 cos(k • ai) cos(k • a 3 )e^ 



r 



(10) 



where the second component turns to be real. The nor- malized factors are given by 
I 



G' n (k) 



1 



^2&(k)4 k + (46(k) 



3e «k) cos 2 (k • a 3 ) + 12c(k)e„ k 
I 



(11) 



The new eigenvectors |w^ k | 

havior at ko; the new singular points brought force by 
are k = ±(-|,^^). Thus according to the two 



recovers the well-defined be- sion of the BZ. Because |it„ k ) and \u' k ) are ill-defined 

only at ko and k , respectively, we are free to deform 
this division as long as ko(k )^V(V'). This corresponds 
to the gauge degree of freedom 2 ^. At keVflV, the two 
choices of wave functions are different by a phase factor 
Kk>=e l9 " (k) Kk), i.e., A / „(k)=A„(k)-V k9n (k), where 



different gauge choices, the BZ is now divided into two 
regions V and V' as shown in Fig. 2. The wave functions 
\u n k) are used onto the region V, while \u' n]g ) apply to 
V'. Note that there remains some freedom in the divi- 



J 



G' n (k) er4 [e Kk cos(k • ai) + 2 cos(k ■ a 2 ) cos(k ■ a 3 )e^] 
Gn(k) i[4 k -4cos 2 (k.a 2 )] 

I 



(12) 



Thus one obtains the values of nonzero Chern number 



for lower and upper bands as follows 

c n = -^- i [A„(k)-A;(k)]-dk : 



— sgn(sin <f) (lower band) 
sgn(sin ip) (upper band) 



I dg n (k) 
2t JdV 

(13) 
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FIG. 2: (Color online). Division of the Brillouin zone of the 
kagome lattice model into two regions V (red area) and V' 
(blue area). 
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which is confirmed by the explicit calculation based 
on the complex-contour integration technique. The 
nonzero Chern number implies that the wave function 
cannot be written as a single function for the entire BZ. 
This also affects the definition of field operators a„k 25 - 
Let a„k denote the annihilation operators when kgV. 
Then \u' nk )=e l9n ^ \u n k) implies the field operator to 
be e _I9 "( k )a n k at keV'. This phase mismatch between 
patches in the BZ provides the quantized Hall conduc- 
tivity. 



III. ORBITAL MAGNETIZATION AND ITS 
EFFECTS 

Now we turn to study the OM and its various effects. 
In the present two-dimensional case, the magnetization 
and Berry curvature become pseudoscalar quantities M c , 
Mq, and fl n . From Eq. (3) one can see that for band 
insulator there will be a discontinuity in OM if the in- 
tegral of the Berry curvature over the entire BZ, or the 
anomalous Hall conductivity, is nonzero and quantized. 
The size of the discontinuity is given by the quantized 
anomalous Hall conductivity times E g /e, where E g is the 
energy gap. Because an analytic derivation of the OM for 
2D kagome lattice is very tedious, thus we turn instead to 
a full numerical representation in terms of the eigenstates 
|u„k) and eigenenergies e n k given in Eqs. (5)-(6). 

We first show in Figs. 3 structures of the energy 
bands e n k, the crystal orbital moment m„(k), and the 
Berry curvature 0„ (k) along high-symmetry lines in the 
BZ for the values of 0=0 (upper panels) and 0=7r/2 
(lower panels). In the case of 0=0, the time-reversal 
symmetry is preserved. From Fig. 3(a) one can see 
that for 0=0 the middle band and upper band touch 
at the corner points (denoted by kx) in the BZ, around 
which the upper and middle bands exhibit a cusp, 63k 
(e2k)=l±V / 3|k — kic|. In this case, the low-energy quasi- 
particle excitations can be well understood within a 
(2+l)-D Dirac fermion field theory. The lower band 
becomes dispersionless (£ik=— 2) at 0=0, which reflects 
the fact that the 2D kagome lattice is a line graph of 
the honeycomb structure^. This flat band touches at 
the r point of the BZ with the middle band, whose 
dispersion around T looks like an isotropic paraboloid, 



FIG. 3: (Color online). (From left to right) band strcture, 
Berry curvature f2„(k), and conventional orbital magnetic 
moment m n (k). The blue, black, and red curves correspond 
to the lower, middel, and upper bands respectively. The spin 
chirality parameter is chosen to be 0=0 for upper panels while 
0=7r/2 for lower panels. 

ff2k=— 2+|fc 2 . Note that the energy spectrum for 0=7r 
(not shown here) is a particle-hole conjugate of that for 
0=0; therefore, the upper band becomes flat with an 
eigenvalue of 2. The Berry curvature fi n (k) and the in- 
trinsic orbital magnetic moment m„(k) for three bands 
are zero everywhere in the whole BZ in the time-reversal 
symmetric case of 0=0, as shown in Figs. 3(b)-(c). This 
can be simply seen by the fact that the wave function 
I link) for 0=0 is real. Generally they satisfy the prop- 
erty m n (— k)=— m n (k), fi n (— k)=— Q n (k) under time- 
reversal symmetry and m n {— k)=m n (k), fi n (— k)=f2„(k) 
under spatial inversion symmetry. Thus in a crystal with 
both time-reversal symmetry and spatial inversion sym- 
metry, fi n (k) and m„(k) will disappear in the BZ. Fig- 
ures 3(d)-(f) show e n k, fi„(k), and m ra (k) for the value of 
0=tt/2. In this case the time-reversal symmetry is bro- 
ken. One can see that there does no longer exist degen- 
eracy in the energy spectrum and two gaps among three 
bands now open. The middle band turns to become flat 
for 0=7r/2, and the energy spectra shows a particle- hole 
symmetry. Note that generally the energy spectra has no 
particle- hole symmetry except for the cases of 0=±7r/2. 
This particle- hole symmetry at 0=7r/2 is also reflected 
in f2 n (k) [Fig. 3(e)] through the feature that the Berry 
curvature for the flat middle band is zero everywhere in 
the BZ, while the Berry curvatures for lower and upper 
bands are nonzero and different with a sign. The intrin- 
sic magnetic moment m n (k) for 0=7t/2 is shown in Fig. 
3(f). One can see that unlike the Berry curvature, m-2(k) 
for the middle flat band is nonzero in the BZ and its am- 
plitude is the sum of Wi(k) and 7713 (k), which are equal 
in the presence of particle-hole symmetry. 

Figure 4(a) shows the OM (JA) as a function of the 
electron chemical potential /i for the value of 0=7r/3, in 
which case neither time-reversal symmetry nor particle- 
hole symmetry is preserved. One can see that initially 
the OM rapidly decreases as the filling of the lower band 
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FIG. 4: (Color online), (a) Orbital magnetization of the 
kagome lattice, and (b) its two components M c (red curve) 
and Mn (blue curve) as a function of the electron chemical po- 
tential /i for (j>=Tv/3. The shaded areas correspond to the three 
groups of bands. To suppress the divergence at band/gap con- 
tacts, we have used the temperature of ksT— 0.05. 

increases, arriving at a minimum at ^=—2, a value cor- 
responding to the top of the lower band. Then, as the 
chemical potential continues to vary in the gap [region I 
in Fig. 4(a)] between the lower band and middle band, 
the OM goes up and increases as a linear function of /i. 
This linear relationship in the insulating region can be 
understood by Eq. (3), from which one obtains 

dM e f d 2 k 

^ = ^E/(^(k) (14) 

occu 
n 

Thus when the chemical potential varies in the gap be- 
tween £ik and £2k, only the lower band £ik is occupied 
and dM/dpi = -{e/h)C x . For 0=tt/3, C x = -1. Thus 
dM/dn = e/h, i.e., the OM linearly increases with the 
chemical potential in the insulating region I as shown in 
Fig. 4(a). When the chemical potential touches the bot- 
tom of the middle band, then the OM suddenly switches 
down and rapidly decreases again when the chemical po- 
tential goes through the middle band. The turning be- 
havior at the band/gap contacts becomes divergent at 
ksT=0. This discontinuity is due to the singular behav- 
ior of O n (k) and m„(k) at the BZ edge points k=k , 
which will play their role when the fc-integral is over the 
entire BZ. When the chemical potential lies in the gap 
between the middle band and upper band, then the OM 
goes up again and increases linearly with /i as shown in 
the insulating region II in Fig. 4(a). Since the Chern 
number of middle band is zero, thus from Eq. (14) and 
Fig. 4(a) one can see that the slope of the OM curve in 
the insulating region II in Fig. 4(a) is the same as that 
in the insulating region I. 

The totally different behavior of the OM in the metal- 
lic and insulating regions, as shown in Fig. 4(a), is due to 
the different roles M c and Mq play in these two regions. 
For further illustration, we show in Fig. 4(b) M c (red 
curve) and Mq (blue curve) as a function of the chemical 
potential, their sum gives M. in Fig. 4(a). One can see 




FIG. 5: (Color online). Orbital magnetization (black curve) 
of the kagome lattice, and its two components M c (red curve) 
and Mn (blue curve) as a function of the spin chirality param- 
eter <j) for the value of \x — — 2, corresponding to the situation 
that the lower band is partially occupied while the other two 
bands are empty. 



that overall M c and Mq have opposite contributions to 
M,, which implies that these two parts carry opposite- 
circulating currents. In each insulating area the conven- 
tional term M c keeps a constant, which is due to the fact 
that the upper limit of the fc-integral of m n (k) is invari- 
ant as the chemical potential varies in the gap. In the 
metallic region, however, since the occupied states varies 
with the chemical potential, thus M c also varies with /x, 
resulting in a decreasing slope shown in Fig. 4(b). The 
Berry phase term Mq also displays different behavior be- 
tween insulating and metallic regions. In the insulating 
region, Mq linearly increases with /z, as is expected from 
Eq. (3). In the metallic region, however, this term sen- 
sitively depends on the topological property of the band 
in which the chemical potential is located. For the lower 
and upper bands with nonzero Chen number, one can 
see from Fig. 4(b) that Mq remains invariant, while for 
the middle band of zero Chern number, it increases with 
the chemical potential fj,. On the whole the comparison 
between Fig. 4(a) and Fig. 4(b) shows that the metallic 
behavior of M. is dominated by its conventional term M c , 
while in the insulating regime Mq plays a main role in 
determining the behavior of M. 

Figure 5 shows the OM M. (black curve) and its two 
components M c (red curve) and Mq (blue curve) as a 
function of spin chirality parameter <f> for the value of 
fj,=— 2. This value of (i ensures that the lower band is 
partially occupied in the whole range of <f>. Thus Figure 
5 describes the metallic behavior of M.. Two prominent 
features can be observed from Fig. 5: (i) The OM is anti- 
symmetric with respect to (f>, M(—<f))——Ai(4>)- This im- 
plies opposite circulating currents carried by Bloch states 
with <p and those with — 0; (ii) In the metallic situation, 
one can see that the amplitude of Mq is much smaller 
than that of M c . As a consequence, the OM is domi- 
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nated by its conventional part M c in the whole range of 
4>. In the insulating region (not shown in Fig. 5), how- 
ever, the amplitude of Mq largely increases and can be 
even larger than the conventional contribution for some 
values of <j> [see Fig. 4]. 

For transport studies of the OM, here we analyze the 
properties of ANE in the 2D kagome lattice. The rela- 
tion between the OM and ANE has been recently found 2 . 
To discuss the transport measurement, it is important to 
discount the contribution from the magnetization cur- 
rent, a point which has attracted much discussion in the 
past. Cooper et ali^I have argued that the magnetiza- 
tion current cannot be measured by conventional trans- 
port experiments. Xiao et al.— have adopted this point 
and built up a remarkable picture that the conventional 
orbital magnetic moment M c does not contribute to the 
transport curent, while the Berry phase term in Eq. (2) 
directly enters and therefore modifies the intrinsic trans- 



port Hall current equation as follows 

JH= - ^ E >< E/ 702/n( r > k ) fi «(k)-V X Mfl(r), 

In the case of uniform temperature and chemical poten- 
tial, obviously, the second term is zero and the Hall effect 
of 2D kagome lattice is featured by nonzero Chern num- 
ber as discussed by Ohgushi et al.— as well as in this pa- 
per. In the following, however, we turn to study another 
situation, where the current-driving force is not provided 
by the electric field (E=0). Instead, it is provided by a 
statistical force, i.e., the gradient of temperature T. In 
this case, Eqs. (15) and (2) give the expression of in- 
trinsic thermoelectric Hall current as j x = a xy (—\7 y T), 
where the anomalous Nernst conductivity a xy is given 
by^ 



1 e x - 



d 2 k 
(2nf 



(e nk - M) fn + k B Tln (l + e -0C**-*O 



(16) 



r 




FIG. 6: The intrinsic anomalous Nernst conductivity a xy 
(divided by the temperature T) of the kagome lattice as a 
function of the electron chemical potential fi for <j)—n/3 and 
fc B T=0.05. The shaded areas correspond to the three groups 
of bands. 



Figure 6 shows a xy of the 2D kagome lattice as a func- 
tion of the chemical potential for 4>—tt/?> and ksT—O.Qh. 
One can see that the ANE disappears in the insulating re- 
gions, and when scanning the chemical potential through 
the bands, there will appear peaks and valleys. Remark- 
ably, a similar peak-valley structure was also found by 
the recent first-principles calculations in CuCr2Se4_ ;c Br a ; 
compound—. The ANE of this compound was recently 



measured by Lee et al£ as a function of Br doping x 
which is used to change the chemical potential /i. Due 
to the scarce data available, until now the peak-valley 
structure of a xy revealed in Fig. 6 and in Refj^ has not 
been found in experiment, and more direct experimental 
results are needed for quantitative comparison with the 
theoretical results. Interestingly, the expression for a xy 
can be simplied at low temperature as the Mott relation 2 , 



klT 



<- xy 



(17) 



where a' xy (fio) is the derivative of Hall conductivity 
with respect to the zero-temperature chemical potential 
(Fermi energy) fiQ. Thus one can see that unlike AHE, 
ANE is given by the Fermi-surface contribution of the 
band structure and Berry curvature. Another unique fea- 
ture of a xy is its linear dependence of temperature. 

Finally, let us consider the response of the OM to the 
external magnetic field B. It should be noticed that the 
above semiclassical theory is carried out up to the first 
order in the external perturbation 4 . Thus our discussion 
is valid in the weak magnetic field. From Eq. (2) one 
obtains the magnetic susceptibility of the OM, 



Xo = 



dM 
~d~B 

E 



(18) 



-x Q c) 



d 2 k 

A o 



m 2 (k) 



dfn 



:Q n (k)m n (k)/ TI 



-(c) 

Obviously, the first conventional term Xo describes the 
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FIG. 7: (Color online), (a) Orbital magnetic susceptibility \o 
of the kagome lattice, and (b) its two components contributed 
from M c (red curve) and Mn (blue curve) as a function of the 
electron chemical potential for cp=iv/3 and feflT=0.05. 

Fermi-surface contribution to \o and thus will disappear 
in the insulating region. Whereas the second term Xo^ 
denotes the Berry phase correction. Clearly, to calcu- 
late Xq 1 ^ all the occupied states within the Fermi level 
should be taken into account. Therefore, the Berry phase 
correction will give a fundamental contribution to Xo in 
the insulating region for the systems with nonzero Chern 
number. 

Figures 7(a)- (b) show the orbital magnetic susceptibil- 
ity and its two components as a function of the chem- 
ical potential, respectively. From Fig. 7(b) one can 

see that the conventional part Xo always displays the 
paramagnetic property, while the Berry phase correc- 
tion contributes a diamagnetic response to B. Again, 
in the metallic regions, the total magnetic susceptibility 
is dominated by its conventional part and thus displays 
the paramagnetic feature. Whereas in the insulating re- 
gions the nonzero Xo comes solely from the Berry phase 
correction and displays the diamagnetic feature. 

Up to now, we have concentrated on the ferromag- 
net represented by the double-exchange model. As Oh- 
gushi et al>22, have pointed out, the present theory is 
also applicable to the ferromagnet based on the Hubbard 
model. Furthermore, the spin-orbit coupling gives the 



spin anisotropies, which introduces the tilting of the spins 
from the perfect ferromagnetic alignment as assumed in 
Eq. (4). The present 2D kagome lattice may be rele- 
vant to the recent experiments on pyrochlore compounds 
R 2 Mo 2 7 (R=Nd, Sm, Gd)iLi£, which are itinerant fer- 
romagnets on the verge of a Mott transition on the py- 
rochlore lattice. It can be expected that easy-axis spin 
anisotropy in these compounds produces the spin chiral- 
ity by the symmetry consideration. Thus the present 
results of the OM properties may be verified in such sys- 
tems. 



IV. CONCLUSION 

In summary, we have theoretically studied the prop- 
erties of OM and its effects in the 2D spin-chiral ferro- 
magnetic kagome lattice. The spin chirality parameter 
<f> in the lattice produces nonzero Chern number and re- 
sults in profound effects on the OM properties. We have 
found that the two parts M c and Mq in OM opposite 
each other, and yield the paramagnetic and diamagnetic 
responses respectively. In particular, due to its Chern- 
number property, the magnetic susceptibility of Mq re- 
mains to be a nonzero constant when the Fermi energy 
is located in the energy gap of the kagome lattice. It has 
been further shown that the OM displays fully different 
behavior in the metallic and insulating regions, because 
of different roles M c and Mq play in these two regions. 
The anomalous Nernst conductivity has also been stud- 
ied, which displays a peak-valley structure as a function 
of the electron chemical potential. The experiments on 
ferromagnets are urgently expected to realize these inter- 
esting theoretical results. 
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